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THE OUTBURST OF THE BLAZAR S5 0716+71 IN 2011 OCTOBER: SHOCK IN A HELICAL JET 
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Abstract 

We present the results of optical (R band) photometric and polarimetric monitoring and Very Long 
Baseline Array (VLBA) imaging of the blazar S5 0716+714 along with Fermi 7-ray data during a 
multi- waveband outburst in 2011 October. We analyze total and polarized intensity images of the 
blazar obtained with the VLBA at 43 GHz during and after the outburst. Monotonic rotation of the 
linear polarization vector at a rate of > 50° per night coincided with a sharp maximum in 7-ray and 
optical flux. At the same time, within the uncertainties, a new superluminal knot appeared, with an 
apparent speed of 21 ± 2 c. The general multi- frequency behavior of the outburst can be explained 
within the framework of a shock wave propagating along a helical path in the blazar's jet. 
Subject headings: galaxies: active — BL Lacertae objects: individual (S5 0716+71) — galaxies: jets 
— polarization 
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1. INTRODUCTION 

The blazar S5 0716+71 is one of most intensively 
studied BL Lac objects. The redshift determinations 
(z « 0.3) are based on measurements of the size of 
a marginally detected host galaxy and a location near 
three galaxies with red shifts of 0.26 (Bvchkov a et al.l 
120061 : iNilsson et all 120081 ). Optical variability has been 
studied by many teams; the blazar exhibits persistent 
activity on both long and shor t (intr a day) timescales, 
as sho w n by, e.g., iWagner et al.l dl996l) ; [Chisellini et al 



19971 ) . IBagen-Thorn et al.l (|2006D and iStalin et al 



20061 120091) and by the resu l ts of several WEBT 



campaigns (iRaiteri et"aTI 120031: lOstorero et"aT1 120061 : 



iVillata et aTl l2008: Che n et al.ll2008j ). Violent polarimet- 
ric variabi lity of S5 07 1 6+71 , among other b l azars , was 
studied bv llkeiiri etail (|2011[ ). lUemura et all (|2010[ ). us- 
ing a Bayesian approach, suggest that the behavior of 
polarization in S5 0716+714 implies the presence of a 
number of polarization components haying q uite short 
timescales of variations. lLarionov et al.1 (|2008h reported 
that a huge optical outburst of S 50716+71 in 2008 April 
was accompanied by a ~ 360° rotation of the position an- 
gle of the electric vector (EVPA) of linear polarization. 

In this paper, we report the results of our observations 
of S5 0716+71 during a major optical outburst in 2011 
October, which coincided wi th an unprecedented 7-ray 
outburst ( Blinov et al.l l20lTh . We present our observa- 
tions in § [2j describe the phenomenological model in § [31 
and apply it to the data and draw conclusions in § [4] 

2. OBSERVATIONS AND RESULTS 

The observations reported here were collected as a part 
of a long-term multi- wavelength study of blazars that fol- 
lows nonthermal activity at optical and 7-ray frequencies . 
An overview of this program is given in lMarscherl {2012). 
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2.1. Optical Observations 

We obtained observational data at the 70-cm AZT-8 
reflector of the Crimean Astrophysical Observatory, 40- 
cm LX-200 telescope in St. Petersburg, both equipped 
with identical photometers-polarimeters based on ST-7 
CCDs, and 1.8-m Perkins telescope of Lowell Observa- 
tory using the PRISM camera. The photometric mea- 
surements in R band were supplemented by observations 
at the 2-m Liverpool Telescope at La Palma, Canary Is- 
lands, Spain. Polarimetry at AZT-8 and the Perkins tele- 
scope was made in Cousins R band, and in white light 
at LX-200, with effective wavelength close to R. The 
photometric errors do not exceed 0^02. 

Instrumental polarization was found via stars located 
near the object under the assumption that their radiation 
is unpolarized. The Galactic latitude of S5 0716+714 is 
28° and Ay = 0^1, so that interstellar polarization (ISP) 
in this direction is less than 0.5%. To correct for the ISP, 
the mean relative Stokes parameters of nearby stars were 
subtracted from the relative Stokes parameters of the 
object. This accounts for the instrumental polarization 
as well. The errors in degree of polarization are less than 
1% (in most cases less than 0.5%), while the EVPA is 
determined with a precision of 1-2 degrees. 

2.1.1. Results 

Our regular photometric and polarimetric monitoring 
of S5 0716+71 started in 2005. Since then, we have 
obtained more than 1000 data points almost uniformly 
spread over a 7-year interval, although we intensified our 
observations during the periods of highest activity of the 
source. 

Figure [1] presents the flux and polarization behavior 
of S5 0716+71 for 2005-2011. We supplement this plot 
with a panel showing the 7-ray light curve from Fermi 
LAT open access data in order to show that the most 
prominent 7-ray activity ever recorded for this source 
was observed during the October 2011 optical outburst. 
In Fig. [5] we show a blowup of the last six months of 2011. 
Table 1 lists the photometric and polarimetric data of 
S5 0716+71 around the time of the outburst. 
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TABLE 1 

Photometry and polarimetry of S5 0716+71 during 2011 October 

outburst 



JD-2400000 R aR p 

(days) (mag) (mag) (%) 



ap EVPA aEVPA Telescope 



(%) 



(°) 


(°) 




246.7 


0.5 


LX-200 


262.9 


0.5 


AZT-8+ST7 


258.7 


0.1 


AZT-8+ST7 


257.6 


2.1 


AZT-8+ST7 


235.6 


0.1 


AZT-8+ST7 


243.0 


1.4 


AZT-8+ST7 


254.0 


1.5 


AZT-8+ST7 


269.8 


1.0 


LX-200 


266.9 


0.4 


AZT-8+ST7 


258.1 


1.2 


LX-200 



55823.516 13.496 

55826.543 13.002 
55827.554 13.012 
55828.547 13.582 

55829.544 13.542 
55830.560 13.598 
55831.574 13.219 
55832.531 13.373 
55833.533 13.361 
55836.488 13.384 



0.012 


20.11 


0.41 


0.003 


8.64 


0.17 


0.002 


16.90 


0.24 


0.003 


5.51 


0.40 


0.004 


16.75 


0.28 


0.003 


9.77 


0.50 


0.002 


3.17 


0.16 


0.003 


11.66 


0.41 


0.003 


20.14 


0.32 


0.003 


19.93 


87 



Note. — Table 1 is published in its entirety in the electronic edition of the 
Journal. A portion is shown here for guidance regarding its form and content. 
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Fig. 1. — From top to bottom: evolution of 7-ray flux density 
(purple points mark upper limits), i?-band flux, and degree and 
position angle of optical polarization of S5 0716+71 during 2005- 
2011. Red horizontal line in the bottom panel corresponds to the 
position angle of the radio jet in 2011. 
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Fig. 2. — Blowup of Fig [T] centered around 2011 outburst. Red 
arrows on R flux panel mark positions of quasi-periodic optical 
flares. 



The optical light curve of S5 0716+71 displays violent 
variability on both long (months-years) and short (days- 
weeks) timescales. Changes in the degree of polarization 
seem to occur erratically: within the time interval 2005- 
2009 it varied from to 15%, while later the mean level of 



polarization increased, reaching a record value of 36% on 
the night of 2011 May 14 (JD 2455696). The EVPA also 
exhibited prominent changes. To solve the ±180° ambi- 
guity, we have added/subtracted 180° each time that the 
subsequent value of EVPA is > 90° less/more then the 
preceding one. 



Outburst of S5 0716+71 in 2011 October. 
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2.2. Radio Observations 

The BL Lac object S5 0716+71 is monitored by the 
Boston University group with the Very Long Baseline 
Array (VLB A) at 43 GHz within a sample of bright 7-ray 
blazar^] with an angular resolution of ~0.1 mas, which 
corresponds to 0.45 pc at the source distance (H = 70, 
A m = 0.3, An = 0.7). The VLBA data were cali- 
brated and imaged in the same manner as discussed in 
Uorstad et all (|2005t ). 

We fit the VLBA data at 16 epochs from 2010 De- 
cember to 2012 May by a model consisting of a number 
of components with circular Gaussian brightness distri- 
butions. During this period we identify 6 components, 
Al, A2, A3, Kl, K2, and K3, in addition to the core, AO 
(Fig. [3J, SJ), which is, presumably, a stationary feature 
located at the southern edge of the jet in the images. 
Figure [5] shows the evolution of the distance of knots 
from the core, while Figure [6] plots the trajectories of the 
knots. 




MJD 55820 55850 55897 55953 55991 56019 56073 



Fig. 3. — Uniformly weighted VLBA images at 43 GHz of moving 
emission knots in the S5 0716+071 jet. Total intensity contours are 
in factors of 2 starting at 0.25%. The restoring beam is identical 
to that shown in Figure [4] New supcrluminal knot K3 appeared on 
MJD 55850, which coincides within 1-cr uncertainty with the time 
of the highest peak in the 7-ray light curve. 



Knots Kl and K3 move to the northeast along 
slightly different trajectories with apparent speeds of 
18.1+0.5 c and 21+2 c, respectively. An apparent 
speed of ~20 c seems to be c ommon in the parsec scal e 
jet of S5 0716+71 (see, e.g., iRastorgueva et all 1201 ih . 
Knot Kl undergoes both an acceleration along the jet 
^11=0.23+0.06 mas yr~ 2 and perpendicular to the jet 
fi± =-0.19+0.06 mas yr~ 2 . An extrapolation of the mo- 
tion of K3 back to the core suggests that it passed 
through the core on MJD 55850+10. Figure [4] shows 
total and polarized intensity images of the quasar while 
K3 emerged from the core. 

Features Al, A2, A3, and K2 are rather stable in dis- 
tance with respect to the core, although they show scat- 
ter in posit ion angle that e xceed s the uncertainties of the 
modeling. Uorstad et al.l (|2001[) found that "stationary 
hot spots" are a common characteristic of compact jets, 
with the majority of such features located within a range 
of projected distances of 1-3 pc from the core. These au- 
thors propose three categories of models for stationary 
components in supersonic jets: a) standing recollimation 
shocks caused by imbalances between the pressure inter- 
nal and external to the jet; b) sites of maximum Dopplcr 
beaming where the jet points most closely to the line 
of sight; and c) stationary oblique shocks, where the jet 

4 http://www.bu.edu/blazars 



bends abruptly. We conclude that knots Al, A2, and 
A3 fall most likely in category a, since they are quasi- 
stationary with an observed "lifetime" of several months 
to a ye ar. In addition, acc ording to numerical simula- 
tions bv I Gomez et al.l Jl997), when a moving knot passes 
through a standing recollimation shock, the components 
blend into a single feature and split up after the collision 
with no changes in proper motion of the moving knot. 
This scenario seems to be observed in the case of Kl 
approaching A2 (Fig. [5J. Knot A3, however, could be 
a trailing compone nt that forms behind a strong shock 
(jAgudo et al.ll200i1 ). since it is observed just after the ap- 
pearance of the fast moving knot K3. The latter might 
change the internal pressure in the jet, resulting in the 
disappearance of knots Al and A2 (see Fig. [5]). Knot K2 
may belong to category c, since the position angle of K2 
is different from that of Al, A2, and A3 (Fig. [6]) and the 
kinematics of moving knot Kl changes in the vicinity 
of K2. In the case of a stationary feature produced at 
a bend in the jet, the proper motion of a moving knot 
approaching the stationary feature should decrease be- 
fore and increase aft er the knot passes through the bend 
(|Alberdi et al.l IT993I ). Such a behavior can be inferred 
from Figure [5] when knot Kl approaches and recedes 
from knot K2. 

2.3. Gamma-ray Observations 

The 7-ray emission from S5 0716+71 was detected 
throughout the entire period of Fermi observations start- 
ing in 2008 August, although most of the time at a 
level < 5 x 10 ph • cm _2 s _1 . We analyze the Fermi 
LAT data for the period of 2011 January-November 
with 2-day binning, as well as with 1-day binning and 
a 6-hour shift, during the huge 7-ray outburst of 2011 
October, when the flux from the source was close to 
2 x 10 -6 ph-cm s . This prevents missing any possible 
short-lived events, to make the correlation analysis more 
robust, and to avoid the dependence of the results on the 
start of the time bins. We use the standard Fermi anal- 
ysis software package Science Tools v9r27pl, with the 
instrument response function P7S0URCE_V6, the Galac- 
tic diffuse emission model gal_2yearp7v6_v0, and the 
isotropic background model iso_p7v6source. The back- 
ground models include all sources from the 2FGL within 
15°. The spectrum of S5 0716+71 as well as spectra of 
all background sources are modelled as a power-law with 
photon index fixed to the catalog values. 

We calculate the discr ete correlation function (DCF) 
ijEdelson fc KrolikJ I1988D of the optical and 7-ray flux 
variations during 2011. The results are given in Fig. [7J 
Optical variations lag those at 7-ray energies by ~ 1.4 
days (DCF centroid position). 

3. MODELING OF PHOTOMETRIC AND 
POLARIMETRIC BEHAVIOR 

We suppose, following Mars cher et al.l (|2008lL that 
most of optical photometric and polarimctric variabil- 
ity arises when a compact emission region (e.g., a shock 
wave) propagates downstream from the black hole, fol- 
lowing a spiral path. (Alternatively, the jet could have 
a helical geometry.) We notice that prominent optical 
outbursts are quite often both preceded and followed by 
minor flares. We interpret this phenomenon as a mani- 
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Fig. 4. — Total (contours) and polarized (color scale) intensity maps of S5 0716+714 at 43 GHz with a total intensity peak of 2.51 Jy/beam 
and polarized intensity peak of 65 mjy/beam; contours represent 0.15, 0.30,...,76.8% of the peak of intensity; the color bar (right) gives 
the color code for the polarized intensity from 5 mjy/beam to 65 mjy/beam. The elliptical gaussian beam is 0.24x0.15 mas 2 (FWHM) at 
PA=-10° (plotted in the bottom left corner); yellow line segments indicate direction of polarization; red and blue solid lines correspond to 
approximate positions of the core AO and superluminal knots Kl and K3. 





0.5 0.4 0.3 0.2 0.1 -0.1 -0.2 
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Fig. 6. — Trajectories of the jet knots. 



festation of the oscillating Doppler beaming of the emis- 
sion ( lighthouse effect, see iCamenzind fc Krockcnbcrgcr 
Ip9a) ). The observed series of outbursts correspond to 
the time intervals when the viewing angle of the shock 



Lag, days 

Fig. 7. — Discrete correlation function between optical and 7-ray 
flux variations of S5 0716+714. Positive lag means that optical 
variations follow those at 7-ray energies. 



wave is at a minimum. 

The main parameters that determine the visible be- 
havior of the outburst are: 

• jet viewing angle 9 

• bulk Lorentz factor T of the shocked plasma 

• temporal evolution of the outburst F' = Fo • 

{^(l|(!lJ)|/2-),*f>to; the factOT k is responsible 
for different time scales of the rise and decline of 
the outburst; primed quantities refer to the plasma 
frame 

• Doppler time contraction in the observer's frame, 
At bs = S 1 • Ai src 



Outburst of S5 0716+71 in 2011 October. 
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• shocked plasma compression rj = ratio of post- 
shock to pre-shock density 

• spectral index of the emitting plasma a 

• pitch angle C of the spiral motion and helical field 

• period of the shock's spiral revolution in the ob- 
server's frame P bs- 

The viewing angle of the compact emission region ip is 
obtained from the relation: 

costp — cos 8 cos C + sin#sin£cos(27rf b s /P bs)- (1) 

Relativistic aberration leads to a change in direction of 
the normal to the shock front ip: 

ip = arctan[sin^/(r(cos<y9- \/l - r~ 2 ))]. (2) 

If the unshocked magnetic field is chaotic upstream of 
the shock wave, the fractio nal polarization of the shocked 
plasma radiation (see, e.g.. lHughes. Aller. fc Allenll985| ) 
is: 

_ a + 1 (1 -rr 2 )sin 2 y> 
P ~ a + 5/3 2 - (1 - 77- 2 )sin 2 V' 

and the position angle is determined by the direction of 
the projected minor axis of the shock wave: 

8 = arctan[Csin£ obs /(Ccos£ obs - 9)}. (4) 

The Doppler factor 8 is given by 

5 = l/(T(l-/3 cos </>)). (5) 

The periods of rotation in the source frame P src and 
observer's frame P bs are connected through the relation 

Pob s = P src (l + z)(l-/3cos0cosC), (6) 

where z is the source's redshift. 

Knowledge of P b s (from the timescale of recurrence of 
outbursts), estimates of (3 and 9, obtained from analysis 
of radio images, and evaluation of the spiral's radius r 
allow us to constrain the values of P src and C and to 
calculate the Doppler-contracted time in the observer's 
frame, i bs- Doppler boosting of the radiation from the 
shocked plasma, as a discrete source, leads to 

Pshock = P s 'hock^ 3+Q - ( 7 ) 

We notice that in most cases when the source is in 
a quiescent state and no outburst is detected, there re- 
mains polarization at the level of < 10%. We assume that 
this is the polarized radiation of the ambient jet, presum- 
ably changing on substantially longer timescales (except 
possibly for random fluctuations due to turbulence in 
the jet). In the local frame the ratio PI = P s 'hock/^j'et 
is of the order of a few percent. The radiation from the 
undisturbed jet is also subject to Doppler boosting, and 
we (rather arbitrarily) set Tj Ct ~ 0.7I\ Since the am- 
bient jet is a roughly steady and continuous source, its 
radiation is amplified by S 2+a . Notice that, in spite of 
the low value of PI, in the observer's frame the boosted 
radiation of the shock may far exceed that of the ambient 
jet. 



We add emission from the quasi-stationary features — 
probably standing shocks — whose fluxes become signif- 
icant after the moving shock crosses them. We postulate 
that the temporal behavior of this added radiation is of 
the same form as that of the shock: 

Apr no p f cxp(-|(t-t +tdci)|/r), t<t 

^jct - rcz- r - ^ exp( _| (4 _ to+tdel) | /feT)) t>tQ 

Here P2 accounts for the ratio of the ambient jet to 
shock luminosity, with time delay tdci- In our model we 
suppose that this radiation is added to the "constant" 
radiation of the jet, but not to its polarized part. Fi- 
nally, the observed radiation of the jet is Pj e t = (Pj' et + 

AF- et )8f+ a and total observed flux P tota i = Pshock + Pjet- 
We consider a single variable source plus an initially 
constant source whose polarization vector is displaced 
relative to the jet direction by an angle Ax (see, e.g., 
lMarsche"rll2010[ ). Accounting for the contribution of this 
"constant" polarized source (F- et and Pj e t)> we get the 
absolute Stokes parameters: 

Q = PshockP C0S(26) + P/ ot <5j t > jo t CO S (2A X ) 

U = PhockP sin(20) + P/ et 4t >jct sin(2A X ) . 
Finally, the normalized Stokes parameters are given by 

q = Ql Ptotai / Q \ 
U = [//Ptotai- W 

Figure [5] shows the dramatically different behavior of 
flux and polarization with all of the above mentioned pa- 
rameters fixed, allowing only the Lorentz factor of the 
shock to change. We see that the model of a shock 
propagating in the helical path may in a natural way 
explain the successions of outbursts, often observed in 
blazar light curves, as a manifestation of a single event. 
A variety of observed patterns of photometric and po- 
larimetric behavior may be explained by the interplay of 
the shock Lorentz factor, jet viewing angle, helical path 
pitch angle, and shape and length of the outburst. We 
do not take into account the light travel delay caused by 
the finite size of the shock wave; if included, this would 
result in some smoothing of the model light curve and 
polarization behavior. 

4. DISCUSSION AND CONCLUSIONS 

Quasi-periodic flares in the light curves of blazars are 
expected in models of magnetically launched jets, where 
the plasma moves with rela tivistic velocities following he- 
lical paths (fvTahakis 2006). In light of this, it is rather 
surprising that numerous extensive observational cam- 
paigns have uncovered only a few cases that could be 
interpreted as a manifestation of a lighthouse effect from 
spiral relativistic motion, e.g., 3C 273 (P = 13 d , opti- 
cal V band ICourvoisier et al.lH988l ); OJ 287 (P = 40 d , 
optical VRI bands I Wu et"all[200l : BL Lac (P = 308 d , 
optic al photometry and polarimetrv IHagen-Thorn et al. | 
2002); and AO 0235+16 (P = 17 d 'in X-ravs lRani et all 
2009). Radio variability tends to occur on timescales 
longer than these periods. However, radio mapping pro- 
vides some picturesque examples of periodic structures 
that correspond to longer timescales in AGN jets, as in 
the case of PKS 0637-752 (2 x 10 3 years iGodfrev et all 
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Fig. 8. — Model evolution of flux and polarization parameters 
with 3 different values of T. 



12012ft . As noted in lLarionov et all (|2010[) . the jet might 
exhibit a fractal helical structure, where geometric scales 
that increase with distance from the central engine cor- 
respond to increasing variability timescal es. Recent po- 
larim etric studies of the nucleus of M87 (Pcrlman et al. 
12011ft have revealed strong variations of the polarization 
vector, synchronous with the flux variations, that can 
be interpreted as a manifestation of a shock propagating 
through a helical jet. 

We have obtained densely sampled optical photomet- 
ric and polarimetric data during an unprecedented 7-ray 
outburst in S5 0716+71, and analyzed the changes in the 
structure of the blazar at 43 GHz. We find a 180° ro- 
tation of the position angle of linear polarization, which 
occurred at the epoch of maximum optical brightness. 
This event would have been missed if our observations 
did not occur consecutively from 3 sites (St. Petersburg, 
Crimea, Arizona) . This implies that such events may be 
rather common during blazar flares. 

Figure |H] shows the time evolution of the observational 
parameters (7-ray and optical flux, degree of polarization 
and positional angle) within the model discussed above. 
Values of the adjustable parameters are listed in TableO 
Several of the trends of the data resemble similar features 
in the model flux and polarization curves. Particularly, 
we see a succession of quasi-periodic flares before and 
after the main outburst (marked with vertical arrows in 
Fig [2]) and fast oscillation of the degree of polarization 
around the main flare. This conforms with the model 
presented here. Turbulence may play a role in modulat- 
ing the behavior (e.g., the discrepant flux and degree of 
polarization near MJD 55675) beyond the capabilities of 
the simplified model proposed here. 

The 7-ray outburst started r?4 before the optical, 
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which corresponds to ~ 0.05 pc distance in the source 
frame. The emergence of knot K3 at the same time as the 
optical/7-ray outburst leads to the conclusion that these 
events were co-spatial. However, we note that our model 
uses a lower Lorentz factor of the shock wave (r = 10) 
and a wider viewing angle {6 = 5.8°) than that obtained 
from the radio images (T = 21). This could reflect ac- 
celeration of the plasma flow in the jet downstream from 
the regio ns where the optica l/7-ray outburst took place 
see, e.g., iHoman et afll2T)09l ). 
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TABLE 2 

Fitting parameters for the photometric and polarimetric behavior of S5 0716+71 in 2011 

October- November. 



e° 


c° 


Pjet,% 


k 


r 


r 


Pohs 


-Psrc 


T 


To 


Rl 


R2 


a 


V 


*del 




(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 
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